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These lectures introduce some of the basic theory of cosmic microwave background 
(CMB) polarization with the primary aim of developing the theory of CMB polar- 
ization from inflationary gravitational waves, as well as some of the related theory of 
weak gravitational lensing (cosmic shear) of CMB polarization. We begin with produc- 
tion of polarization by Thomson scattering. We then discuss tensor-harmonic analysis 
(the "grad-curl" or "E-B" decomposition) on a flat and full sky in some detail. The 
Boltzmann/Einstein equations required to predict the CMB temperature/polarization 
pattern due to primordial gravitational waves are derived. We show that gravitational 
waves produce a curl component of the CMB polarization while density perturbations 
(at linear order) do not. We then show how cosmic shear induces a curl component from 
a curl-free surface of last scattering. We describe, though in less detail, how higher-order 
correlations can be used to subtract the cosmic-shear-induced curl. Several exercises are 
provided. 

1. Introduction 

Just five years ago, the small-scale structure of the CMB was hidden behind a 
veil of experimental limitation. That is no longer the case. We now know empiri- 
cally that the CMB has a wealth of detailed information in its temperature [1] and 
polarization pattern [2]. We are now also confident in our theoretical understand- 
ing of these fluctuations: they are produced by a nearly scale-invariant spectrum 
of primordial density perturbations. With a solid theoretical foundation, we are 
poised to move considerably further with future more precise CMB experiments. 
Prospects for new advances include an even more detailed picture of the physics 
at the surface of last scattering, the physics of inflation, gravitational-wave back- 
grounds, the distribution of mass in the intcrmcdiate-redshift and current Universe, 
and perhaps the behavior of the dark energy. 
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Fig. 1. Radiation incident on an electron from the ±x directions scatters into the z direction by 
shaking the electron in the ±y direction. 



These lectures will not discuss all these wonderful possibilities. Instead, we 
focus here on detailing some of the basics of the theory of CMB polarization with 
the aim primarily of understanding the polarization due to inflationary gravitational 
waves, as well as the related effects of weak gravitational lensing ("cosmic shear") 
on CMB polarization. There is now a vast literature on CMB theory, including 
a number of reviews [3] and even a few books [4]. Thus, experts are unlikely to 
find anything new in here. However, some students may find that the discussion 
and level of detail for the specific subjects on which we focus provide a useful 
complement to these other more comprehensive sources. In particular, we discuss 
tensor harmonics (the "curl-grad" or "E-B" decomposition) on both a flat and 
a full sky in detail, and we provide a detailed justification, from scratch, of the 
statement that gravitational waves produce a curl component in the CMB while 
density perturbations do not. We show in detail that cosmic shear induces a curl 
component in the CMB polarization, but our discussion of how this shear may be 
subtracted with higher-order correlations is a bit sketchier. We do not deal with 
any of the many important issues in data analysis. 

The plan of the lectures is as follows: We begin in Section 2 by reviewing how 
Thomson scattering produces polarization. Section 3 develops tensor harmonic 
analysis on a flat sky and Section 4 deals briefly with flat-sky correlation functions. 
Section 5 deals with tensor harmonic analysis on the full sky and Section 6 with 
the correlation functions. Section 7 provides some comments, and in Section 8 we 
discuss some of the effects of cosmic shear on the polarization pattern. A few exer- 
cises are sprinkled throughout to help conscientious students develop their facility 
with subject matter. We reference primarily the most important original papers on 
relevant subjects and/or those that were particularly useful to us in the preparation 
of these lectures, but our referencing does not provide a full guide to the literature. 
For that we refer the reader to the more comprehensive reviews [3] . 

2. Polarization of CMB: Theory 

Let us begin by understanding in the most basic terms why the CMB should be 
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Fig. 2. Radiation incident on the electron with equal intensities in the ±x and ±y directions will 
be unpolarized when scattered to the z direction. 



polarized [5] . Consider scattering of an electromagnetic wave propagating in the ±£ 
direction in the plane of the page that then scatters to the observer (perpendicular 
to the page) by 90°; the scatterer is a free electron (Fig. 1). If the incident radiation 
is linearly polarized with an E vector in the z direction, then the electron oscillates 
in the z direction and emits no dipole radiation toward the observer. But if the 
incident radiation is polarized in the y direction, the electron oscillates in the y 
direction and the dipole radiation the observer sees is linearly polarized in the ±y 
direction. If the incident radiation is unpolarized then only the component with E\\y 
is scattered toward the observer, and the resulting scattered radiation is polarized in 
the y direction. Likewise if the incident radiation comes only from the ±y direction, 
then the observed scattered radiation will be polarized in the ±x direction. 

Now if we have radiation incident from the ±x and ±y directions with equal 
intensities (see Fig. 2), then the polarizations in the ±x and ±y directions will 
cancel, and the scattered radiation will be unpolarized. But if the intensity in the 
x direction slightly exceeds that in the y direction, then the scattered radiation will 
be slightly polarized in the y direction. 

To proceed further we must now introduce the Stokes parameters Q and U for 
linear polarization. Most generally, a monochromatic electromagnetic wave propa- 
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Fig. 3. The linear-polarization states described various combinations of Stokes parameters Q and 
U. 

gating in the z direction has an electric-field vector, 

E x = a x cos(wt - £ x ); E y = a y cos(u>t - £ y ). (1) 
The Stokes parameters are then the intensity, 

I = al + a 2 y , (2) 

linear-polarization parameters, 

Q = a 2 x -a y , U = 2a x a y cos(^ - (3) 

and circular-polarization parameter, 

V = 2a x a y sm(£ x - £ y ). (4) 

The last Stokes parameter V vanishes as Thomson scattering induces no circular 
polarization (but see Ref. [8]). The parameter Q quantifies the polarization in the 
x-y directions while U quantifies it along axes rotated by 45°. Fig. 3 illustrates the 
linear polarization described by various combinations of Q and U. 

If we rotate the x-y axes by an angle a about the line of sight z, then the new 
x'-y' coordinates are 



x \ _ I cos a sin a \ / x 
y' J \ — sin a cos a J \ y 



but the Stokes parameters (Q, U) transform as [9] 

Q' \ _ f cos 2a sin 2a \ f Q 
U' / { -sin2a cos2a j { U 



(5) 



(6) 



Formally, (Q, U) are two quantities that under a coordinate transformation, 

x' i = A\x k , (7) 
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transform as 

P[ 3 =A\Ayp kl . (8) 
More explicitly, (Q, U) are the components of a symmetric trace- free 2x2 tensor, 

Q U \ _^ ( cos a sin a \ ( Q U \ ( cos a —sin a ^ . . 

(9) 



U —Q J V —sin a cos a J \ U —Q J \ sin a cos a 

or alternatively and equivalcntly, a spin-2 field. 

We now return to the generation of polarization by Thomson scattering. Gener- 
alizing our earlier heuristic results, the magnitude of the polarization of the scattered 
radiation is proportional to the magnitude of the quadrupole of the radiation inci- 
dent on the scattering electron, and the orientation of the polarization is determined 
by the orientation of the quadrupole. More precisely, from the angular dependence 
[6], 

of Thomson scattering, where <tt is the Thomson cross section and e' and e are the 
polarization of the incident and scattered radiation, respectively, it can be shown 
that the polarization of radiation scattered from an electron cloud of optical depth 
t«1 into the z direction is 



Q- l U=^ra 22 . (11) 

Here, a 2 2 is the radiation quadrupole moment incident on the electron cloud. More 
precisely, a 2 i is the coefficient of the spherical harmonic Yii{Q,§) in a spherical- 
harmonic expansion of the incident-radiation intensity in a coordinate system in 
which the line of sight is the z direction, x and z are in the scattering plane, and Q 
and U are measured with respect to the x and y axes. (See Ref. [7] for expressions 
for more general orientations.) 

Exercise 1. Verify Eq. (11). 

Finally (!) we can see why the CMB must be polarized. We see angular vari- 
ations in the temperature of the CMB implying temperature (and thus intensity) 
inhomogeneities at the CMB surface of last scattering. Thus, most generally, a 
given scatterer at the last-scattering surface will see an anisotropic distribution of 
incident radiation leading to polarized scattered radiation. For example, the radia- 
tion in the center of Fig. 4 will be polarized in the direction of the two cold spots 
in the temperature field. 

3. Harmonic analysis for Q and U on a flat sky 

According to inflation (as well as other structure- formation theories), the tem- 
perature pattern on the sky is a single realization of a random field. The simplest 
(and for a Gaussian field only non-trivial) statistic to describe the temperature 
pattern is the power spectrum C, TT . 
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Fig. 4. The radiation scattered from the center of this quadrupolar temperature pattern will be 
polarized as shown. 

Suppose that we now have in addition to the temperature pattern the polar- 
ization, quantified by Stokes parameters Q(n) and U(h) measured as a function of 
position h = (9, (f>) on the sky. Again, since Q and U are components of a symmetric 
trace-free (STF) 2x2 tensor, they depend on the coordinate system we choose. To 
discuss physics, we will thus want to find a coordinate-system-independent repre- 
sentation of this tensor field. Later, we will need to do this on the full sky, but as a 
warmup, we will first do the simpler case of a flat sky (which also serves as a good 
approximation to a small region of the sky). 

Before moving further, it should be noted that the formalism we are about to 
develop also applies to cosmic-shear (weak gravitational lensing) maps, where the 
cllipticity parameters e + and e x are, like Q and U, components of a STF 2x2 
tensor [10]. 

Once we have measured the polarization, Q{9) and U (9), as a function of position 
9 = (6 X , 6y ) on a flat region of sky, we have measured the polarization tensor field, 



which is sometimes written as a vector, P a (Q,U) = P(Q,U), although this is a 
dangerous thing to do as Q and U do not transform as the components of a vector. 
The polarization is also sometimes written as a complex number, 




(12) 



P= \P\e 



(Q 2 + U 2 ) 1/2 e 2 



Q + iU, 



(13) 
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where | P\ = (Q 2 +U 2 ) 1 ^ 2 is the magnitude of the polarization and a = (1/2) arctan([//Q) 
is its orientation relative to the x axis. 

We now define gradient ('G', also known as E) and curl ('C, also known as B) 
components of the tensor field that are independent of the orientation of the x-y 
axes as follows: 

V 2 P G = d a d b P ab ; V 2 P C = e ac d b d c P abl (14) 

where e ab is the antisymmetric tensor. 

We can write more explicit expressions for these G and C components in Fourier 
space. Writing 

P ab (0) = J ^2 p a4)e- lt \ (15) 

p ab (i) = / d 2 9P ab (ey r -°, (16) 

the Fourier components of Pq (9) and Pq (0) are 

5 l {ll-ll)Q{l) + 2l x l y U{l) 

m) = 2 WTTi ' (17) 

Pc{l) - 2 zfTTf • (18) 

Exercise 2. Verify these expressions for the Fourier components P G and Pq, and 
verify that they are invariant under a rotation of the 9 X — 6 y axes. 

For a temperature map T{9), the temperature power spectrum C ; TT is defined 
from 

lf{l)f{V))^{2n) 2 5{l + V)Cr, (19) 

where the angle brackets denote an average over all realizations. Likewise, the 
statistics of the polarization field are determined by polarization power spectra 
Cf G , C cc , and Cp c defined by 

(P G {l)P G {f)) = (27r) 2 S(l + l ) )C^ G , (20) 
(pc(f)Pc(?)) = (2tt) 2 S(1 + 1 ) )CP c , (21) 
(p G (f)P c (F)\ = (2k) 2 5{1 + T')C? c . (22) 



If we also consider cross-correlation of the polarization with temperature, then there 
are in total six power spectra, 

■X^X^V)) = (27r) 2 d(l + V)Cf' x \ (23) 

where X u X 2 = {T, P G , P c }. 
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9i 

Fig. 5. Correlation functions for polarization are best denned for Stokes parameters Q r and U r 
measured with respect to axes aligned with the line connecting the two points being correlated. 

Now suppose we have a given temperature/polarization map and then consider 
a parity inversion; e.g., consider a reflection about the x-axis. Then 

Oy ► 6y, Q > U ► £/, l x > Ixi (24) 

Also, 

f(l)^f(f), P G (f)^P G (f), P c (l) -> -P c (l). (25) 

In other words, G and T are parity even, while C is parity odd. Thus, if the physics 
that gives rise to T/P fluctuations is parity conserving, we expect 

C; TC = C ; GC = 0, parity invariance (26) 

in which case the statistics of the T/P map is determined entirely by the four power 
spectra, CJ T , C ; TG , Cf G , and C cc . 

4. Correlation functions 

For a temperature map, the correlation between the temperature T(9i) and 
T(9 2 ) at any two points is 

(T$)T(&)) = C TT (Vi - 2 |) , (27) 

and the correlation function Ctt(#) depends only on the distance 9 between the 
two points. The same is not true for Q and U . Since Q and U are not rotational in- 
variants, their correlations will depend also on the orientation of the line connecting 
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the two points. Instead, one can define rotationally invariant correlation functions 
by considering the components Q r and U r of the polarization defined with respect 
to the line connecting the two points (Fig. 5). We can then define the correlation 
functions, 

(Qr(0l)Qr(O 2 )) = C QQ (|&-&|), (28) 

(£/ r (0i)*7 r (0 2 )) = Cuu (|& - &|) , (29) 

(<2 r (0i)*7 r (0 2 )) = C QU (|&-&|). (3°) 

Since Q — > Q an d U ^ —U under a parity inversion, Cqu = if parity is con- 
served. If we correlate the polarization with temperature, then there are another 
three (TT,TQ,TU) correlation functions, one of which (TU) vanishes if parity is 
conserved. 

The correlation functions can be written in terms of the power spectra, 

f°° 1 HI 

C QQ {9) + Cuu((9) = -/ — [Cf G + Cf c ]J (W) ) (31) 

Jo n 

fOO 1 TT 

c QQ (o) - c vv {9) = - _[q GG -q cc ]j 4 (Zfl), (32) 

Jo ^ 

where J v (x) are Bessel functions. These relations can also be inverted to give 
power spectra in terms of correlation functions. For derivations and details, see 
Refs. [10,11]. 

Finally, the correlation functions at zero lag give us the mean-square polariza- 
tion, 

<P 2 ) = (Q 2 + U 2 ) = J™ l -^[Cr + CF C ] - + (P 2 c). (33) 

Exercise 3. Write an expression for the temperature power spectrum Cj T in terms of 
the temperature autocorrelation function C TT (0). Now do the same for the gradient 
and curl power spectra in terms of the QQ and UU autocorrelation functions. 



5. Harmonic analysis on the full sky 

If our maps extend beyond a small region of the sky, we will have to come to 
terms with the fact that the sky is a curved surface. Moreover, as we will see below, 
theoretical predictions for the power spectra will require solution of Boltzmann 
equations that require that we keep track of the distribution function for photons 
moving in all directions. 

We thus generalize the tensor Fourier analysis (G and C) that we did above for 
STF 2x2 tensors that live on the 2-sphcrc. Our discussion in this Section follows 
Rcf. [11]; a different but equivalent formalism is presented in Rcf. [12]. In the usual 
spherical polar coordinates 9, (f>, the sphere has a metric, 

sin%)- ^ 
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The polarization tensor P a t must be symmetric P a b = Pba and trace-free g ab P a b = 0, 
from which it follows that, 

rab(Tl) - 2 ^ -U{n)sm9 -Q(n)sin 2 9 J ' 

(The alert reader will note that the sign of U has changed here relative to Eq. 
(12). This is simply because for the full sky we are following the slightly different 
conventions used in Rcf. [11].) The factors of sm.9 also follow from the fact that 
the coordinate basis (9, 4>) is orthogonal but not orthonormal. 

Let us now review some of the rules of differential geometry on the two-sphere. 
Covariant derivatives of scalar, vector, and tensor fields are 

S-. a = S^ n = r, 6 + y c C T ab :c = T a \ c + T db T a cd + T ad T b cdl (36) 

and ':' denotes covariant derivative, 5 a = (dS/dx a ), and the Christoffel symbols 
are ^ 

rtc = ^9 ad (9db,c + 9dc.b - 9bc.d)- (37) 

Also, 

o'-ab yj2t-72 q , ridb q i ^ TD-d q yy2 q — n:a r> — r?c r> — r?a 

<J afc— VVD+it D-.db +2 ' = a = acb: = °" 

(38) 

Since the sphere has no boundary, we can integrate by parts, 

j> d 2 n^X ab Y :ab = -j> d 2 h^X ab :a Y :b = j> d 2 h^X ab :ba Y, (39) 

where Y(n) is a scalar function and § d 2 h denotes integration over the sphere. The 
antisymmetric tensor is 

Cat = ( J ) > ( 4 °) 

and 

£cae C fc = g a b = —tac^b, tabbed = 9ac9bd — 9ad9bc, £ab:c = 0. (41) 

Also, e ab T a b = for a symmetric tensor, T a b = Pba- 
For the two-sphere, the Ricmann tensor is 

11 1 

Rabcd = ^ R tabbed, Rab = -^Rg a b, £° Rabcd = R^cd, e"' Rabcd = ^R^bd- 

(42) 

For two STF tensors M ab , N ab (i.e. g ab M ab = g ab N ab = e ab M ab = e ab N ab = 0), the 
second of Eqs. (41) gives us 

M ab N cd e ac e bd = -M ab N ab . (43) 

Also, 

g = \ 9ab \=^e, et=(_l e S[ ; 9 ), (44) 
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T% = - sin 9 cos 9, T+ H = Y% = cot 9, (45) 
and all other components vanish. Then, for a scalar function Y(h), 

Y-gg = Y :0 g, 

Y..0^ = Yoj, - cot OYrf, 

Y-M = Y i<t><t> + sin9 cos 9Y, e . (46) 
A symmetric rank-2 tensor M ab has 'divergence', 

M ab .. ab = M 6e . ee + 2M e \ H + M++ M - sin0cos0M^, e 

+ 2 cot 8M ee . e + 4 cot 0M**^ + (1-3 cos 2 9)M <t " t ' - M ee , (47) 

and 'curl', 

M ab .. ac e c b = sin9(M^. eg + M^^e)-cscd(M 0e . g<p + M^.^)-cot9csceM e0 ^ 
+5 cos 9M H fi + 3 cos 9M^ A + 3 (cos 9 cot 9 - sin 9) M 6tt> , (48) 

Exercise 4. Verify Eqs. (47) and (48). 

Any STF 2x2 tensor field can be written as the 'gradient' of some scalar field 
A(n), ' 

A-.ab - ^9abA' c C , (49) 

plus the 'curl' of some other scalar field B(h), 

l(B :ac e c b + B :bc e c a ). (50) 

Just for comparison, a vector field is analogously decomposed as 

V a = V a A + e ab V b B. (51) 

Since any scalar field on the sphere can be expanded in spherical harmonics (e.g. 
for the temperature), 

f = i+EE *U Y ™w> (52) 

l— 1 m— — l 

where 

4rn) = Y o J dnT(n)Y{ lm) (h), (53) 

it follows that the polarization tensor can be expanded in terms of basis functions 
that are gradients and curls of spherical harmonics, 

V ab (h) ' ' 



EE [4rn) Y (f m )ab(n) + af lm) Y^ m)ab (n)] . (54) 



T 

1=2 m=-l 
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The expansion coefficients are given by 

4m) = Y Q j dflV ^{n)Y^*{h), af lm) = ±rj dnV ab (h)Y^*(n), 

(55) 

and 

Y (L)ab = N l (^(Imy.ab ~ ^9abY( hn):c C ^j , (56) 
Y (?m)ab = -y (yY(lm):act C b + Y(l m):hc € C , (57) 

constitute a complete orthonormal set of basis functions for the G and C components 
of the polarization. The quantity, 



is a normalization factor chosen so that 



/ dnY^ )ab {n) = J dnY^ )ab (n) Y^h) = 6 w 5 mm >, (59) 

JdnY { f r : )ab (h)Y { f,^ ) (n) = 0. (60) 
Also, we can integrate by parts to write alternatively, 



«(L) = y q I dhY( lm) {h)V ab - ab {h), (61) 



at r 

4m) = yJ dfl Y ^ {fl) Vab ' ac ^ e "- ^ 

Finally, since {T, Q, U} e 3?, 

4m) = {- l ) m 4,-mV (63) 

where X = {T,G,C}. 

Exercise 5. Consider a vector Geld V a that lives on the surface of a two-dimensional 
sphere. Show that this vector field can be written as 

oo I 
1=1 va— — l 

where Y^ m ^ a and Y^ m ^ a are gradient and curl vector spherical harmonics. Derive 
expressions for these harmonics and show that they are orthonormal. 

The temperature/polarization power spectra are now 

4m)4'm')) = hl'dmm', (64) 



Cabella and Kamionkowski 13 




pure 

G 



pure 

c 



Fig. 6. Examples of G (top) and C (bottom) polarization patterns. 

for X, X' = {T, G, C}. Of these six, two (TC and GC) are zero if parity is conserved. § 
For future reference, the Yn m ) ab and Y^ m ^ ab are explicitly, 



Y {lm)ab {h) = \ . . , (65) 

2 V X (lm) (n)sm9 -W {lm) (n) sm 2 8 J 



and 



where 



q ,^_Nif W { i m) {n) X {lm) {h)sm.( 

lm)ab\ n ) — ~Y 

c iV, / -X(im)(fi) W {lm) (n) sin t 
Y {lm)ab {n) = — , (66) 

2 V W(im)( n ) Sln ^ A" (im) (n) SU1 Z 6> 



/a 2 a m 2 \ / <9 2 \ 

^(M(») = [q02 ~ cote Q- e + —2^) Y (im)(n) = (^2— + i(/ + 1)J K (Im) (n), 

(67) 

and 

07,777, i O \ 

X{im){n) = -— e {a- e - cotd ) (68) 
In terms of the spin-2 harmonics ±2^(;m) used in Rcfs. [12,13], 



W (Im) (n) ± iX (lm) (n) = W |±|jj ±2 F (;m) . (69) 



§Our GG,CC, and TG moments can be identified with the EE, BB, and TE moments of Rcf. 
[12] through a£> = af m /V2 and a^ m j = a ^ m l^- Also, the (7j here reduce in the small-angle 
(large-Z) limit with those in Section 4 as long as the angles in the flat-sky limit are given in radians; 
it is left as another exercise to the reader to verify this statement. 
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Note that if we replace (Q, U) by (U, -Q), then G -» C and C -» -G. This tells 
us therefore, that a pure-G polarization pattern becomes a pure-C pattern if we 
rotate each polarization vector by 45°, and vice versa. Examples of G and C type 
polarization patterns are shown in Fig. 6. 

6. Correlation functions on the sphere 

Correlation functions on the two-sphere are defined analogously to those on 
the flat sky. The line connecting the two points is now replaced by the great arc 
connecting the two points [11]. To calculate the two-point correlation functions, it 
is convenient to choose the north pole of the coordinate system to be one of the 
points. The temperature autocorrelation function is then 



For polarization, we need a bit more algebra. We begin by deriving expressions for 
the functions that appear in our definitions above of the tensor harmonics: 





Iml'm' 




(70) 




(71) 




(72) 



with 





(74) 



The Legendre functions have the following asymptotic limits, 



If 1 (cos 0) 



(_l)(m+|m|)/2 (; + | m |)l 

2\ m \\m\\ (l-\m\)\ 




9 







(m ^ 0) 



(75) 




(76) 
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We also note that at the north pole, Xi m (9 = 0) 7^ and Wi m {6 = 0) ^ only for 
to = 2: 



wwo,o) = l f l t 1( {i-Z {Sm ' 2 + Sm '- 2) ' 



and 



i 2/ + l(H_2)! 
A (;m) (0,0) = 2 V n _ 2)! ( to2 ~ ^ ' 



The QQ autocorrelation function is then defined as 
C?QQ(g) / Qr{ni) Q r (n 2 y 



-*0 -M3 / ft 1 .f l2=C os0 



(77) 



(78) 



(79) 



where Q r is measured with respect to the great arc connecting the two points. Let 
Q(0, 0) = lim ^o Q(0,0). Then, 



or 



C QQ (#) 



cQQ{9) = ^ Q(o,o)cmo) ^ 



Qr(0,0) Q r (0,O) 



(80) 



Iml'm 



]T NtNvl [af lm) W {lm) (0,0) - af lm) X {lm) (0,0) 

I'm' \ 

[ a fl*m') W *l'm')(8>°) - a fl* m ') X (l'm')(°^) 



= E 



21 + 1 

8tt 



JVi[Cf-(W(i2) + ^(1-2)) + iCf-(X ( * I2) - X$,_2))]> 



07 I 1 

= E -^^ 2 ( cos + cfGr 2 (cos ^)] , 



2tt 



where we have used in the first line, 



(81) 



Q(n) - 2P M (n)=T ^ ^ ^^^(nj-og^ln)], (82) 

Z=2 m=-Z 

C/(n) = -2 esc 6V <p(h) 

OO I 

= " T E E N l \ a (lm)X(lrn) (n) + fl^ m) W(, m) (ft) 



(=2 m=-( 

Similarly, the UU correlation function is 

2Z + 1 



C VV {6) = ]T [CfGUcosO) + C?G n (cos9)] , 



(83) 



(84) 
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and the TQ correlation function is 

C^(0) = (^pl9lM) =£^cr G Pf(c 0S S). (85) 



T o To / ni . A2=cos9 4tt 

Note again that the TU and QU correlation functions vanish if parity is conserved. 

Exercise 6. If parity is broken, do the TU and QU correlation functions vanish? If 
not, calculate them in terms of the (parity-violating) power spectra. 

7. Calculation of Predicted Power Spectra 

Now that we have figured out how to describe the polarization field properly, 
let's proceed to see what theoretical models (e.g. primordial adiabatic density 
perturbations and gravitational waves) predict about the polarization. Since these 
calculations are extremely involved in practice, here we only outline the calculation. 
We will, however, be able to see precisely that gravitational waves produce a curl 
component, while density perturbations do not. 

To begin, we follow the pioneering paper of Polnarcv [14] to derive the angular 
distribution of photon intensities in the presence of a gravitational wave (GW). To 
begin, we suppose that the photons do not scatter. In this case, the photon energies 
are affected only by the form of the metric. Let us consider a single monochromatic 
plane-wave gravitational wave, which appears as a tensor perturbation to the FRW 
metric, 

ds 2 = a 2 {n)drf - a 2 (n)[dx 2 (l + h+) + dy 2 (l - h+) + dz 2 }, (86) 
where -q is the conformal time and 

h+(x,ri) = h{ri)e thn e- tkz , (87) 

describes a plane wave propagating in the z direction. This is a linearly-polarized 
gravitational wave with "+" (rather than "x") polarization. Here h(n) is the am- 
plitude; at early times when krj < 1, h(rj) ~const, but then h(rj) redshifts away 
when k-q > 1 (see, e.g., Ref. [15]). 

Exercise 7. Show that the amplitude h + (rj) solves 

h+ + 2-h + + k 2 h + = 0, 
a 

where the dot denotes derivative with respect to conformal time (actually, anisotropic 
stresses give something on the right-hand side, but we will neglect them). Then solve 
this equation for the case of pure matter domination and for pure radiation domina- 
tion. (Your solution should be consistent with a gravitational-wave energy density 
pew oc a ~ 4 with the scale factor a(t) of the Universe when the mode wavelengths 
are smaller than the horizon.) Think about what the spectrum of gravitational 
waves should be in the Universe today, assuming an initially scale-free spectrum. 
You should find that the slope of the spectrum changes at the scale that entered 
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the horizon at matter-radiation equality. The spectrum is plotted in Fig. 3 of Rcf. 
[16], and you should be able to understand the result plotted there analytically. 

With the metric [g a p = a 2 diag(l, —(1 + h + ), —(1 — h + ), —1)] , the zeroth com- 
ponent Po of the photon four-momentum is the photon energy E multiplied by the 
scale factor P = aE; P is constant if h + = 0. We take h + <C 1 for small-amplitude 
gravitational waves. 

Using the geodesic equation, 

pa= dx^ jl_^d^J_^dr)_d_^ p0 d_ ^ 
dX ' dX dX dx° dX dr\ dr\' 1 
and g^vP^P" = for photons, we arrive at the Sachs- Wolfe effect for this spacetime 
(to first order in h), 

where e x = (1 — /i 2 ) 1 / 2 cos0 and e y = (1 — (J. 2 ) 1 ^ 2 sin (f> are the components of 
the direction of the photon momentum, and /i = cos 9 and <p describe the photon 
direction. 

Exercise 8. Verify Eq. (90). 

Replacing Po by the comoving photon frequency v = Po = av. p h ys =const (in the 
unperturbed FRW spacetime), we find that the gravitational wave induces angular 
intensity variations of the form, 

^ = (91) 

Now let us consider Thomson scattering of these photons by electrons. Scattering 
will change this angular distribution and it may induce polarization as well. To 
include the effect of Thomson scattering, we (following Polnarev 1986 [14]) consider 
an alternative set of "Stokes parameters", and describe the state of the radiation 
propagating in any given direction n by the "vector" , 




1= \ h , (92) 



where Ig — a 2 ,, 1$ = a|, Q = Ig — 1$, and I = 1$ + 1$, and Q and U are measured 
with respect to the 9-(j) axes on the plane tangent to the sky at any given direction 
n={6,4>). 

The distribution function / for photons is now also a "vector" . For example, in 
a homogeneous universe, it is 

/ = /oM ( 1 | =/o(M), (93) 
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where 

/o(") = £h J T I ( 94 ) 

is the usual blackbody distribution function, and the form (1 1 0) indicates no 
polarization. In the presence of the gravitational wave, there will be a perturbation, 
as determined above, so 



/(M) = /o 









/o 


[(:) 











(95) 



Note that /i is a three-component "vector" , and it is most generally a function of 
conformal time r], position x, frequency v, and photon direction (0, (f>) = n. 

As the Universe expands, photons get re-distributed in frequency, polarization, 
and direction by the redshift due to the gravitational wave and also by Thomson 
scattering. Until now we have considered only the gravitational redshift and ne- 
glected scattering, which we now proceed to include. The time evolution of the 
photon distribution function is determined by the equation of radiative transfer, 
essentially the Boltzmann equation for the photon distribution function, 

dr\ drj dx l dv drj 
where g = axn e a is the scattering rate, n e the electron density, and 

where 

o / A* V 2 cos 2(0' - cj>) -i? cos 2(0' - 4>) n 2 n' sin 2(0' - 4>\ 
P=-\ -m' 2 cos 2{4>' - 4>) cos 2{4>' - (f) -it' sin 2(0' - 0) I (98) 



4 



2^' 2 sin 2 (</>'- 4>) 2/j,sm2(<j>' - <j>) 2 W ' cos 2(0' - (f>) 



is the scattering matrix. Eq. (96) says that the total time derivative of the photon 
distribution function (written as an explicit time derivative plus the time evolution 
due to photon motion) is given by gravitational redshift (the first term on the right- 
hand side) in the perturbed spacetime plus the change due to Thomson scattering 
(the second term on the right-hand side) of photons from electrons. The quantity 
J is the angular intensity-polarization distribution that arises after the radiation 
has been Thomson scattered once from an initial distribution /. The scattering 
matrix P looks messy, but it simply tells us how an initial intensity-polarization 
pattern is re- arranged after Thomson scattering once from unpolarized electrons. It 
is straightforward to derive it from the angular-polarization dependence of Thomson 
scattering, 

(x|e-ef; (99) 
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see Refs. [14] and [9] for details. 

The gravitational wave imprints at first an angular intensity pattern, 

ldv = _h 1 _2s Qs2(j)e -ik z d ik 

This intensity pattern is then altered when the photons undergo Thomson scatter- 
ing. Before Thomson scattering, the intensity-polarization distribution function is 
the unperturbed distribution f plus a small perturbation proportional to 

a= i ^ 1 ) (1 -m 2 )cos2</i, (101) 

due to the gravitational wave. After Thomson scattering, the angular intensity- 
polarization distribution is altered; the scattering matrix above then introduces a 
second component to the perturbation proportional to 



1 / (1 + ^ 2 )cos20 
'l + ^ 2 )cos" 
4:fi sin 2(p 



b=-\ -(l + ^ 2 )cos2</> I , (102) 



One might guess that subsequent scatterings would add further more complicated 
angular dependence to the distribution function. However, as first noted by Polnarev 
[14], Thomson scattering of a distribution function of the form b returns a distribu- 
tion function of the form a. In other words, the basis functions a and b form a closed 
basis under Thomson scattering. This is simply a consequence of the time-reversal 
invariance of the scattering process: If Thomson scattering converts a distribution 
a to a distribution 6, then it should turn a distribution b into a distribution a. 

The linearized solution to the Boltzmann equation, including fully the effects 
of gravitational redshift and Thomson scattering, in the presence of a gravitational 
wave must therefore be of the form 

/ = /o + e- 4fcz+4fc "/i, (103) 

where the perturbation is 

/i = a.(r],vo,n)a + (3{r), vo,^)b. (104) 

Here, a(r], Uq, /u) and /3(?7, v$, /x) are coefficients that must be determined by solution 
of the Boltzmann equation. 

Defining an "anisotropy" A = £(1 — y?) with £ = a + /3, and "polarization" 
n = (3(1 + ii 2 ) (which is nonzero only if there is polarization), the radiative-transfer 
equation can be re-written, 



Jo avo v ' 

3 rO\0 ^ 1 



$+[ik(l-iM)+g]0 = ^g I d/ZKl + M 72 ) 2 /?--^!-// 2 ) 2 ]. (105) 
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The first equation generates a temperature fluctuation from the gravitational wave, 

and the second generates polarization through scattering of that temperature fluctuation. 11 

Eqs. (105) are a set of coupled partial differential equations for the time (77) 
and angular (/i) dependence of the distribution functions ^(77, /x) and /3(n, u). In 
practice, these are solved numerically by Legendre transforming £(77, /i) — > (,i(v) 
and (3(ri,n) — > Pi(r)) through 

1 f 1 

6 fa) = g / i d/^fa.^O*), (106) 

&{v) = lj dn/3(ri,n)Pi{n). (107) 

The partial differential equations for £(77, /x) and /3(r/, /j) then become an infinite set" 
of coupled differential equations for £i(rj) and A(^) which are propagated numeri- 
cally from some very early time (where the solutions are analytically determined in 
the so-called tight-coupling approximation) to the present time. 

Exercise 9. Derive the differential equations for £i(r]) and /3;(r?); this is not neces- 
sarily an easy problem. If you want to go further, you can determine the early-time 
solutions to these equations using the tight-coupling approximation, in which the 
scattering rate g — n e OTa is assumed to be huge. 

We will not discuss the numerical techniques (which are highly nontrivial) here, 
but will show results a bit later. For our purposes, we simply need to know that the 
angular polarization pattern induced by this gravitational wave can now be written, 

Q(M) = Jt o ^(2Z + 1)P ; (cos#)(1 + cos 2 0)cos20£, 
l 

{7(0, 0) = ^T ^(2Z + l)P;(cos0)2cos0sin20£;. 

We then get a polarization tensor, 

P o6 (M) = ^£(2Z + l)i5(cosl% 
8 1 

/ (1 + cos 2 0) cos 20 -2 cot 9 sin 20 

X \ -2cotflsin20 -(1 + cos 2 9) esc 2 9 cos 20 

Exercise 10. Verify Eq. (108). This should be easy. 

^Note that the factor ik(l — fi) that appears in Eqs. (105) often appears elsewhere simply as ikfi. 
The reason traces back to the fact that Polnarev writes h-y. oc h(r])e lk ' n and / = /o + f 1 e~ lkz + lkr i 
while other authors usually write/14, oc h(rj) and / = /o + We thank J. Pritchard for clarifying 
this point. 

llThese are truncated at some large value of I; the truncation procedure is not trivial; see, e.g., 
Ref. [17]. 



(108) 




Cabella and Kamionkowski 21 



If we now expand Eq. (109) in tensor spherical harmonics, the resulting tensor- 
harmonic coefficients are 



o? m = \Ni + 1)6 J dnY{ m {n)M?-y. ab {i 



(110) 



which after a little bit of algebra becomes 



Likewise, 



^{6 m2 + 6 mt - 2 )y/2n(2l + l) 

(z+2)(z+i)6_ 2 w+m i(i - m+2 

(21 - 1)(2Z + 1) (2Z + 3)(2Z - 1) (2Z + 3)(2Z + 1) 



= l - Nl ^J + mj dhY^(h)M^ ):ac (h)4 



(111) 



(112) 



2tt 



4 V (2Z + 1) 



_ 2 )[(Z + 2)6-i + (Z-l)6+i]- (113) 



We have thus shown explicitly that both the G and C components are nonzero for 
a gravitational wave. 

We get the contributions to the power spectra Cp G and C cc from this particular 
gravitational wave (in the z direction with '+' polarization) from 



21 + 
16 



+ 1 ^ 1 



G |2 



(Z + 2)(Z + 1)6-2 , 6Z(Z + 1)6 



+ 



+ 



1(1 - 1)6+2 



(2Z-1)(2Z + 1) (2Z + 3)(2Z-1) (2Z + 3)(2Z + 1) 



(114) 



and similarly for C cc . Summing over all Fourier modes, / d 3 k/(2Tr) 3 , and over 
both polarization states, the final result for C GG is 



C, GG 



16tt/ 



fc 2 dfc 



(Z + 2)(Z + l)6- 2 6Z(Z + 1)6 Z(Z - 1)6+2 

(2Z-1)(2Z + 1) (2Z + 3)(2Z-1) (2Z + 3)(2Z + 1) 



(H5) 

and similarly for C GC . Note, finally, that the cross-correlation power spectrum 
vanishes 



C 



GC 



m = l n G* n C 
u lm u lm 

21+1 



m— — l 



o. 



(116) 



as it should. This is because a G m) oc (5„ h2 + S m ,-2), while a^ m ^ cx (8 m .2 - 5 m ,-2)- 
Now what about scalar (density) perturbations? By following steps similar to 
those above, one can show that they do not produce a curl component in the CMB 
polarization. Here we only sketch the calculation and leave out details. Again, 
consider a single Fourier mode of the density field in the z direction. Then the 
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Sachs- Wolfe effect induces an intensity variation proportional to (cos 2 9 — 1/3), so 
now 

/i = a{ri,v ,ii)a + P{r),v ,ij)b, (117) 

with 

V J 

We thus find that for k\\z, U(h) — 0, so 

I 

-l 

sin^ 9 

One finds af m ^ 0, but af m = 0. This follows because Mf^ c e% = which follows since 
M°^ c is diagonal and independent of </>. Therefore, if one detects a curl component 
in the CMB, it is a signature for primordial gravitational waves (see Fig. 7). 

Exercise 11. Show that if the gravitational-wave background is composed entirely 
of gravitational waves with right-handed circular polarization (and no left-handed 
gravitational waves) that Cj G ^ 0; i.e., that there is a cross-correlation between 
the CMB temperature and the curl component of the polarization. (Hint: Consider 
a single circularly-polarized gravitational wave propagating in the z direction; from 
Ref. [18].) 



M$(0,<f>) = sin 2 0P,(cos0) I J ). (119) 



8. Comments on power spectra 

The power spectra shown in Fig. 7 show a variety of features. The origin of 
these features, as well as their dependence on cosmological parameters, has been 
the subject of much study. We will not discuss them in detail but refer the reader 
to the reviews and books listed in the Introduction for more details. Here, we point 
out only a few interesting and relevant features. 

• There are acoustic peaks in the GG power spectrum for density perturbations. 
These are out of phase from the peaks in the TT power spectrum. The 
temperature fluctuation is due primarily to density fluctuations at the surface 
of last scattering, and secondarily to the peculiar velocity at the surface of last 
scattering. The peaks in Cj TT are due to the density perturbations; the troughs 
arc filled in by peculiar velocities, which oscillate (just as in an ordinary 
harmonic oscillator) out of phase by 90°. 

Before recombination, the electrons and photons are tightly coupled, so the 
electrons see no quadrupolc in the photon intensity and thus produce no po- 
larization. Just before recombination, the photons begin to have longer path 
lengths as they begin to decouple. Since this nonequilibrium process depends 
on the time derivative of the baryon density (just like the velocity does), the 
polarization is out of phase with the density, and in phase with the velocity 
[19]. 



Cabella and Kamionkowski 23 




Fig. 7. Power spectra of temperature and polarization (G and C) and cross-spectrum TG for scalar 
perturbations (dotted line) and perturbations due to gravitational waves (solid line). Note that 
the power spectrum of 'curl' component C, due to scalar perturbations, is missing, and the dashed 
line corresponds a reionized model with optical depth r = 0.166 to the surface of last scattering. 
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The power spectra due to gravitational waves drop precipitously for I > 100. 
This is because on smaller scales, the gravitational waves have entered the 
horizon and had time for the their amplitudes to redshift away by the time of 
recombination. 

Reionization induces a large bump in C GG , C GC , and C ; TG at I < 10 [20]. 
This is simply due to scattering of the quadrupole by reionized gas. A rough 
estimate for the I of this peak can be obtained by assuming an Einstein-de- 
Sitter Universe and noting then that 

(z • \ 1/2 1 

Uon~ 200 p^ ^T~ 10 fOT ^cio„~10, (120) 

\ ^rcc / VO 

where I <~ 200 is the acoustic-peak location, and we have used the approxima- 
tion that the sound speed c s ~ c/ \/3 in the primordial plasma is dominated 
by radiation. 

C GG from density perturbations peaks at I <~ 1000, as opposed to I <~ 200 for 
C ; TT : there is more power on small scales in the polarization. This is because 
polarization induced by a particular Fourier mode of the primordial density 
field depends on the gradient of that density field. 

The amplitude of C CG is proportional to the square of the amplitude of the 
gravitational-wave background, which is fixed by the height of the inflaton 
potential during inflation. This can be quantified by 



where V is the inflaton-potential height during inflation, and cJ T ' tensor j s 
the tensor (gravitational-wave) contribution to the temperature quadrupole. 
Since the latter is ~ 10~ 10 , we already known that the energy scale of inflation 
is V 1/>4 < 3 x 10 16 GeV. The curl component of the CMB polarization induced 
by gravitational waves is thus proportional to the scale-height of inflation: 

Cf c oc V. (122) 

Therefore, detection of a curl component in the CMB due to gravitational 
waves would provide not only a "smoking gun" for inflation, but it would also 
tell us the height of the inflaton potential during inflation, and thus provide 
some hints as to the new physics responsible for inflation. In practice, any 
realistic CMB polarization experiment will be able to detect the GW-induced 
CMB curl component only if V > 10 15 GeV; in other words, only if inflation 
took place at the energy scale of grand unification. 

C cc from gravitational waves peaks at I ~ 100 or at <~ 2°. If there is no other 
sources of a curl component (see below) , then detection of the curl component 
is not cosmic- variance limited. 
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• Gravitational waves contribute a roughly scale-invariant temperature power 
spectrum C ; TT at low I, just like density perturbations. However, cosmic 
variance removes our ability to detect an excess of < 10% over the density- 
perturbation contribution to Cj T , even if we could predict theoretically the 
density-perturbation amplitude. 

• If there is no other source of a curl component, and if the large- angle polar- 
ization bump due to rcionization is not big, then the optimal survey strat- 
egy [21] entails a deep integration on a ~ 5° x 5° patch of sky, with an 
angular resolution < 1° fwhm. This is the strategy of experiments like BI- 
CEP and QUIET. We now know from the WMAP detection of a large-angle 
temperature-polarization correlation that the reionization bump is probably 
quite big. In this case, an all-sky experiment with roughly the same instrumen- 
tal sensitivity (noise-equivalent temperature) might have similar sensitivity to 
inflationary gravitational waves. 



Exercise 12. Scale- invariant spectra of density perturbations and of gravitational 
waves each produce a nearly scale-invariant spectrum, 1(1 + 1)C ; ' tcnsor ^constant, 
in the CMB at large scales (i.e., I < 20). Suppose that most of the observed 
anisotropy at large scales comes from density perturbations, and pretend that we 
knew precisely the amplitude of the density-perturbation power spectrum from mod- 
cling smaller-scale fluctuations, other observations, or divine inspiration. If so, there 
would be a limit, set by cosmic variance, to the smallest gravitational-wave ampli- 
tude (quantified by, e.g., 6Cj T,tcnsor , the gravitational-wave contribution to the 
temperature quadrupolc moment) that could be determined from the measured 
CMB power spectrum. Estimate the smallest detectable value of cJ T ' tonsor _ You 
will need to start by showing that the la cosmic-variance error with which each 
Cj T can be measured is ^2/(21 + l)Cf T . 



9. Gravitational Lensing (Cosmic Shear) and CMB Polarization 

Above we showed that density perturbations do not induce a curl in the po- 
larization, and thus concluded that detection of a curl in the CMB polarization 
automatically implies detection of gravitational waves. However, that derivation 
assumed only linear perturbations, in which each Fourier mode of the density field 
is considered independently. If more than one Fourier mode of the density field 
is considered, then different modes can interact and produce a curl, even without 
gravitational waves. Since <5 <~ 1CP 5 at the last-scattering surface, this density- 
perturbation-induced curl component should be small. 

We do know, however, that the observed CMB temperature-polarization map 
will be distorted by cosmic shear (CS), gravitational lensing by large-scale mass 
inhomogeneities between us and the surface of last scattering. The effect of cosmic 
shear is to displace the temperature and polarization from a given direction 8 at 
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the surface of last scattering to an adjacent position, 9 + 59: 

0)= \ Q ] {6 + 56)~ \ Q | (9)+S9-V\ Q j (0). (128) 

where <50 = Vtp is the cosmic-shear displacement angle, and 

V{n) = -2 r dr n M r, hr) (124) 

J dA(r,0)d A {ri s ,0) 

is the projection of the gravitational potential $(r, hr) (obtained from the mass 
distribution through the Poisson equation) along the line of sight, and cU(r i,^) is 
the angular-diameter distance corresponding to the comoving radial coordinate r\ 
by an observer at r^- 
Noting that 

p _ p 21 I 

TT-lf = cos WTH = sin (125) 

'a "r 'j/ l :r ' L y 

our previous transformation [cf., Eq. (18)] between (Q,U) and (G, C) is (using G 
and C as a shorthand for Pq and Pc, respectively), 

G _ 1 ^ cos2(^- sin20j- \ f Q 



C /^2Uin2^ -cosJirA? J«- (126) 

Let us suppose the polarization field at the surface of last scattering has no curl; 
then this relation can be inverted to give 

Q(l) = 2G(f) cos 20 f , U(T) = -2G(f) sin 2<j> T , (127) 

and 

VQ(9) = -2i j G (0 cos20 r r e - irf , (128) 

and similarly for VC/ (9) with cos — ► — sin. The displacement angle is likewise 

500) = -ij -IL V (J) e-* S t (129) 
Thus, the perturbation to Q and U induced by gravitational waves is 

6Q(9) = VQ • V^ = J AL e - itS {VQ ■ V<p) T , 
where 

5Q(l) = (VQ • V^) r = 2 y ^[fi • (f- f!)]G(ri) V (r- fi) cos 2^-, (131) 
<5£/(f) = (VU • Vy>) r = -2 y 70j[^i ' (i"*- fi)]G(fiMf- f^sh^-. (132) 



(130) 
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Although the original map had (by assumption) no curl, the lensed map does: 
C(f) = i[sin20 r Q(f) -cos 20,-17$] 

= J (£^ [fl ' (f ~ ZW&M*- sin2(^— fi ) 

= / (2^ [fl ' (f ~ MM^M^ sin2 <%- ( 133 ) 

Thus, our earlier claim that detection of a curl component constitutes detection of a 
gravitational-wave background is not entirely valid, as we have just shown explicitly 
that a curl component can be induced by cosmic shear [22]. 

Using the power spectrum (|<Pf| 2 ) = Cf v for the projected potential, the power 
spectrum for the curl induced by lensing is, 

c?G = j (£^ [fi • (f - fi)]2 sin2 2 wri^- (i34) 

In Fig. 8 we show numerical results for the CMB polarization curl induced by 
cosmic shear (lensing) [23,24]. If the gravitational- wave amplitude (solid line) is 
large enough, the CS curl will not interfere with detection of gravitational waves, 
and the gravitational-wave signal can be distinguished from the lensing signal by 
the shape of the power spectrum. However, if the gravitational-wave amplitude 
is small (< 4 x 10 15 GeV), lensing (dashed line) produces a background, and the 
gravitational-wave curl cannot be detected [23,24,25]. 

Fortunately, something can be done to help separate the CS-induced curl from 
the GW-induced curl. Seljak and Zaldarriaga [26], Hu and Okamoto [27], and others 
have shown that measurement of higher-order correlations induced by lensing can 
be used to reconstruct 69(6), the displacement, as a function of position on the 
sky. Refs. [23,24] have then evaluated how well these may be used to reduce the 
CS-induced curl. This subtraction is somewhat involved technically, and still under 
active investigation. We give a very brief description here to provide the basic flavor 
of the technique; readers are referred to the original papers for more details. 

In the absence of lensing each Fourier mode of the T (or Q or U) field is statis- 
tically independent: 

(T(T)T(P)} = for J"V ?. (135) 

However, if there is lensing, an observed Fourier mode T(l) has contributions from 
all pairs of temperature and projected-potcntial Fourier modes T(l{) and (fiih) that 
have I = l\ + ?2- Thus, with lensing, 

(T{l)T{?)) = f{l, T)<p{L) j> t (136) 

in the presence of some fixed projected potential <p(0) with Fourier components 
<fi(L). Here, 

f(l,l > ) = Cr(L-l) + Cr(L-P). (137) 
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Fig. 8. Contributions to the CMB polarization power spectra. The upper (blue) solid "gravita- 
tional waves" curve is the curl component due to gravitational waves in the presence of reionization 
with an optical depth r = 0.17 (note the large-angle bump) and the associated (blue) dashed curve 
is that with no reionization. The amplitude of this curve is for the largest inflaton-potential height 
(V ~ 3.5 X 10 16 GeV) allowed by COBE; note that C ; CC ' GW <x V, so the amplitude of this curve 
will be reduced accordingly if V is reduced. The short-dash (green) "lensing" curve is the curl 
power spectrum induced by cosmic shear (weak gravitational lensing due to density perturbations 
between us and the surface of last scattering) . The red "scalar" curve is the GG power spectrum 
due to density perturbations (with reionization), shown here for reference. The dotted "lensing" 
(green) curve is the cosmic-shear contribution to the curl component that comes from structures 
out to a redshift 2 = 1, and the green "lensing" dot-dash curve is the residual cosmic-shear 
power spectrum left after subtraction with higher-order temperature-polarization correlations, as 
described in the text. From Ref. [23]. 
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(There are analogous expressions for polarization, but for simplicity, we deal here 
only with T.) To determine a given L component of tp, we simply sum over all 
pairs T(li), T(l 2 ) with L = l\ + l 2 . In practice we have to weight these pairs 
taking into account noise from the random-field nature of the fluctuations, as well 
as instrument noise. Doing so, Hu-Okamoto find the following estimator for the 
Fourier components of the displacement angle, 

56{L) = / ^T(h)T(l 2 )F(l u l 2 ), (138) 

-l 

, (139) 

and C; 11 ' 1 is the total observed (signal plus noise) power spectrum. Thus, with 
these estimators, the projected potential can be determined as a function of po- 
sition across the sky from the measured temperature map. We can then use this 
projected-potential measurement to reconstruct the polarization pattern at the sur- 
face of last scattering from the (lensed) polarization pattern that is observed. Sim- 
ilar estimators that use the lensed polarization (rather than temperature) can also 
be constructed [27,28], but we do not include them as the expressions rapidly be- 
come unwieldy. The precision with which (p(L) can be reconstructed depends on 
the number of small-scale coherence patches in the temperature-polarization map 
that can be used as 'sources' with which the shear can be reconstructed. Thus, 
high angular resolution and high sensitivity are required. Since the polarization 
power spectrum peaks at I <~ 1000, rather than I <~ 200, there are more small- 
scale coherence patches in the polarization than in the temperature. As a result, 
a high-sensitivity and high-resolution polarization map will be required for optimal 
lensing reconstruction. The degree to which the curl component induced by cosmic 
shear can be reduced depends on how well this reconstruction can be accomplished. 
For plausible assumptions about a post-Planck CMB polarization satellite (such as 
those that NASA is now considering), the CS-induced curl can be reduced by a fac- 
tor <~ 10 in power-spectrum amplitude (as shown in Fig. 8). This requires a full-sky 
map of the temperature and polarization with an angular resolution 0f w hm ^arcmins 
and detector noise-equivalent temperature s ~ 1/iK-^/sec. Higher-order correlations 
may improve upon this subtraction by a bit more [29,30]. 

Exercise 13. Show that the estimators given in Eqs. (138) and (139) are the 
minimum-variance estimators for the projected potential that can be constructed 
from the temperature map. You may need to go to the original papers to see how 
this is done. 

The best strategy to detect inflationary gravitational waves is now a subject of 
some study. If a curl search probes inflaton-potential heights V 1 ^ > 4 x 10 15 GeV, 
then CS-induced shear does not constitute a background. In this case, one strategy 
is to integrate deeply on <~ 5° x 5° patch of sky with moderate (<~ 0.5°) resolution. 
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f(hM) 

2C TT,t c TT,t> 



A(L) 



f(hM)F(hM) 
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For a one-year experiment with NET s > 10/iKVsec (e.g., Planck, BICEP, QUIET), 
this is the best strategy. However, if the large-angle temperature-polarization cross- 
correlation is as large as WMAP indicates, then there will be a large reionization 
bump in the curl power spectrum. If so, then a similar sensitivity to gravitational 
waves might be achievable with a full-sky map with ~ 0.5° angular resolution and 
similar detector sensitivity. 

However, if the effective detector sensitivity is s < lOyuK-^/s, then inflaton- 
potential heights V 1 / 4 < 4 x 10 15 GeV can be probed. For these smaller-amplitude 
gravitational-wave signals, the CS-induced shear will be a foreground for the I ~ 
50 — 100 GW-induced Cp c . This foreground must be removed with higher-order 
correlations, which requires a full-sky map of the temperature and polarization 
with #f w hm ^arcmins. This will have to be the route for post-Planck s ~ l/iK-v/sec 
experiment such as CMBPOL. With quadratic estimators, the lowest detectable 
inflaton-potential height V is V 1/4 ~ 10 15 GeV [23,24], although one may be able 
to do a bit better with higher-order estimators [29,30]. Anything smaller will be 
lost in the CS-induced curl, even after subtraction with higher-order correlations. 
However, with a large-angle reionization bump, as suggested by recent WMAP mea- 
surements, the gravitational-wave signal at large angles may be as easily detectable 
as the smaller-angle signal, without being confused by the CS-induced curl [31]. The 
"best" survey strategy is still not entirely clear, and will depend on experimental 
factors as well as these more theoretical considerations. These questions are the 
subject of several ongoing NASA mission concept studies. 

10. Closing Comments 

In these lectures, we have provided details of some of the theory of CMB polar- 
ization. The reader who has successfully gone through all of the Exercises may take 
pride in the knowledge that he/she has achieved a mastery of the technical aspects 
of the subject comparable to that of many researchers in the field. Of course, these 
lectures still leave many essential aspects of CMB polarization theory uncovered. 
These include techniques for solving the Boltzmann equations and a more intuitive 
understanding of the features of the power spectra. The reader should also be aware 
that the full- or flat-sky tensor-harmonic formalism must necessarily be altered to 
analyze real CMB maps. Such maps will have sky cuts and so the measured regions 
of sky will be less than the full sky and may be irregularly shaped. Sophisticated 
techniques for dealing with a cut sky have now been developed, and many other 
issues that accompany analysis of real maps are highly nontrivial and completely 
neglected here. 

The past few years have been quite exciting for the CMB. There is also clearly a 
very active foreseeable future in the field with a number of targets for experimental- 
ists: e.g., the polarization autocorrelation function, detection of weak-lensing, and 
the longer-term goal of detecting inflationary gravitational- wave background. If the 
past is any guide, however, chances are that the most exciting discoveries are those 
that theorists have not yet anticipated. 
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Appendix: Rotational invariance of an 'spin-2 field' 

Here we demonstrate the rotational invariance of Eq. (14). Let us begin by 
considering the usual rotational matrix, 

4 = A{ Xj . (140) 

The derivative operator transforms as 



OF _ dF dx\ 
dxj dx\ dxj 



From Eq. (140), we obtain 



(141) 



(142 » 



and Eq. (141) becomes 

dx-=dx> Ai - (W3) 
Thus, the derivative operator in the rotated system transforms as 

ir {A ~ 1] -^ <144) 

Taking into account Eq. (8), Eq. (14) becomes 

d[d' 3 P[ 3 = (A-^dkiA-^A^AiP^ = 8 l n 5^d k dtP mn = OkdiPu- (145) 
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